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Abstract 

We perform the threshold resummation for Higgs pair production in the dominant gluon fusion 
channel to next-to-next-to-leading logarithmic (NNLL) accuracy. 

The calculation includes the matching to the next-to-next-to-leading order (NNLO) cross sec¬ 
tion obtained in the heavy top-quark limit, and results in an increase of the inclusive cross section 
up to 7% at the LHC with centre-of-mass energy Ecm = 14 TeV, for the choice of factorization 
and renormalization scales = Q, being Q the invariant mass of the Higgs pair system. 

After the resummation is implemented, we estimate the theoretical uncertainty from the per¬ 
turbative expansion to be reduced to about ±5.5%, plus ~ 10% from hnite top-mass effects. The 
resummed cross section turns out to be rather independent of the value chosen for the central 
factorization and renormalization scales in the usual range {Q/2,Q). 
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1 Introduction 


The experimental collaborations ATLAS and CMS have discovered in 2012 a new particle at 
the LHC 0i. whose properties are so far compatible with the Standard Model (SM) Higgs 
boson. In order to be able to distinguish between the SM and other new physics scenarios, it is 
crucial to measure its properties as accurate as possible. In particular, to understand the relation 
between the discovered boson and the electroweak symmetry breaking mechanism it is essential 
to determine its couplings to fermions and gauge bosons, and verify its proportionality to the 
particle masses. Furthermore, it is necessary to measure the Higgs self-interactions, in order 
to start reconstructing the scalar potential, which is responsible for the spontaneous symmetry 
breaking. While the Higgs quartic coupling is currently out of reach [^, several studies have 
recently shown that a measurement of the Higgs self-coupling can be achieved at a luminosity 
upgraded LHC via Higgs pair production [4 11 . This could be one of the main goals of the 


high-luminosity run of the LHC. 

As it happens for single Higgs production, the SM Higgs pair production is dominated by 
the gluon fusion mechanism mediated by a heavy-quark loop. At leading order (LO) in QCD 
perturbation theory 12 14 , this process can occur either via a box diagram, gg ^ HH, or a. 
triangle diagram, gg —)• H* HH, being the latter the only one sensitive to the Higgs trilinear 
coupling. Given that this process is already one-loop level at LO, higher order corrections are 
very difficult to compute in the full theory. 

The next-to-leading order (NLO) corrections have been computed in the large top-mass ap¬ 
proximation in Ref. 15 , and more recently also the next-to-next-to-leading order (NNLO) cross 
section became available within the same effective theory 16 . The QCD corrections were found 


to be large, resulting in about a 100% increase from LO to NLO, and a still sizeable 20% increment 
from NLO to NNLO at a collider center of mass energy of 14 TeV. The higher order corrections 
are almost completely dominated by soft and virtual terms, as it was shown in Ref. 16 . The theo¬ 


retical uncertainty arising from missing higher orders in the perturbative expansion was estimated 
to be about ±8.5% at NNLO for this energy. On top of this, one should add the uncertainties 
of the strong coupling and parton flux determination, plus the unknown hnite quark-mass effects. 
In order to reduce these uncertainties, one should compute higher order corrections and hnite 
top-mass effects, respectively. 


In this work, we compute the dominant effect of the uncalculated higher-order terms by ex¬ 
ploiting the resummation of soft-gluon emission, working within the large top-mass approximation. 
We provide numerical results for the LHC up to the next-to-next-to-leading logarithmic (NNLL) 
accuracy. Our calculation consistently includes the matching to the NNLO cross section. 

The paper is organized as follows. In Section we dehne our notation and present all the 
expressions needed to perform threshold resummation for Higgs pair production up to NNLL. 
In Section we present the numerical results for the LHC, comparing the hxed order and the 
resummed cross sections in several aspects. Finally, in Section]^ we present our conclusions. 
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2 Threshold resummation for Higgs pair production 


We consider the production of a Higgs boson pair via top-quark loops. We work within the 
large top-mass approximation, where the effective gluon-Higgs coupling is given by the following 
Lagrangian 
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— I Gh - GhH —V 


H 




( 1 ) 


where v ~ 246 GeV is the Higgs vacuum expectation value and stands for the gluonic field 
strength tensor. The perturbative expansions of G// 17-20 and Ghh 
order needed for our calculation, i.e. O (ag). 


21,22 are known up to the 


We start by setting the notation for the fixed order calculation. The hadronic cross section for 
a center-of-mass energy of the collider sh, differential in the Higgs pair invariant mass Q, can be 
written in the following way 


da 


= (J{sh,Q^) I dXi dX2 fa/hAxi,fil) fb/h2{x2,fJ^l) 


a,b 


X / dz 6{ z — 
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X 1 X 2 J 


do ^ Gabiz; tts(/iR), Q Q /i^f) > (2) 


where r = Q'^/sh, and are the renormalization and factorization scales respectively, and do 
is the Born level partonic cross section. The parton densities of the colliding hadrons are denoted 
by fa/h{x, fijp) and the subscripts a,b label the type of massless partons (a, 6 = g,qf,qf, with 
Nf = 5 different flavours of light quarks). The hard coefficient function Gab can be expanded in 
terms of powers of the QCD renormalized coupling as(/i|.) (in the following, the dependence of 
as on /iR is understood): 


Gaj)(x; as, Q^//ip) — (^) • (3) 

n=0 

We use the MS scheme for the renormalization of the strong coupling. 

Given that the soft-gluon resummation has to be carried out in Mellin (or Wmoment) space 
, we introduce now the corresponding notation. We consider the Mellin transform crAr(Q^) 
of the hadronic cross section. The Wmoments with respect to r = sp at fixed Q are dehned 
as follows: ^ 

o-7v(Q^) = [ dr a{sH,Q‘^) ■ (4) 

Jo 

In A^-moment space, Eq. ([^ takes the following simple factorized form 

O'N-liQ'^) = do ^2 fa/hi,N{.g^F) fh/h 2 N{.g^V) Gafe^ Ar(as, Q^//i|’) , (5) 

a,b 

where we have introduced the Wmoments of the parton distributions and the hard coefficient 
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function as 
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Once these A^-moments are known, the physical cross section in z-space can be obtained by Mellin 
inversion by 
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27ri 
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where the constant Gmp that dehnes the integration contour in the At-plane is on the right of 
all the possible singularities of the integrand, as dehned in the Minimal Prescription introduced 
in 
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We want to consider the all-order summation of the enhanced threshold (z —>■ 1) contributions, 
which corresponds to the limit N ^ oo in Mellin space. Given that gg is the only partonic channel 
which is not suppressed in this limit, we only need to consider its contribution. The resummation 
of soft-gluon effects is achieved by organizing the partonic coefficient function in Mellin space as 


— Cgg{as, Q'^/Q'^/fJp) 

■ A^(«s, QV^4li + 0{l/N) , (9) 


The large logarithmic corrections (that appear as Og ln^"^“^ N in Mellin space) are exponentiated 
in the Sudakov radiative factor Ajg, which depends only on the dynamics of soft gluon emission 
from the initial state partons. It can be expanded as 

\nAN(as,\nN; = InA^ g^^\^oas\nN) + g^‘^\Poas in / nl) 

+ as 9^^\f3oas i^N,Q‘^/fi%-,Q‘^/g,l) 

H-oo 

+ 9^^^ (/^oas In N, /ii\, j■ (10) 

n=4 


The function In A resums all the leading logarithmic (LL) contributions ag In""^^ iV, 
contains the next-to-leading logarithmic (NLL) terms ag In"' N, as collects the next-to-next- 
to-leading logarithmic (NNLL) terms ag^^ In" N, and so forth. All the perturbative coefficients 
required to construct the and gi^) functions are known and only depend on the nature of 


the initiating partons. Their explicit expression can be found, for instance, in Refs. 26,27 


The function Ggg{as) contains all the contributions that are constant in the large-A^ limit. 
They are produced by the hard virtual contributions and non-logarithmic soft corrections, and 
can be computed as a power series expansion in as: 


Cggias,Qygl;Q^/gl) = 1 + X] (^) gli QVgl) 

n=l 


( 11 ) 
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The Cgg coefficient, needed to perform N*LL resummation, can be obtained from the N*LO 
hxed order computation. The only process-dependent contribution to Cgg arises from the virtual 
corrections, given that the soft contributions are universal. In fact, in Ref. 28 we derived a 
universal formula for the coefficients needed up to NNLL accuracy, only dependent on the infrared 
regulated one and two-loop amplitudes (more recently the calculation was extended to N^LL in 
). Then, we can obtain the expressions for Cgg^ and C^g from the explicit two-loop 
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ation performed in Ref. 30 . Specifically, we have 
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fin 


O-Q 
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where 'fE is the Euler number, (n represents the Riemann zeta function and (3o = (HC^ — 
2A^/)/127r. The infrared-regulated one and two-loop corrections dg^ and dg^^ are defined by the 
following set of formulae (for fiR = fip = Q) 
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fin 




CTq 


\Clo\ {t+ — t_)Jt_ 
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with t± = — I — 2Mjj =p Q\/Q^ — and where for the sake of brevity we refer the 


expressions 


of 7lP‘\ and to Ref. 30 (supplemented with Ref. 1^ for C^^^ 


HH) 


Finally, in order to proht from the hxed order calculation we implement the corresponding 
matching. As usual, we expand the resummed cross section to 0{ag) subtract the expanded 
result from the resummed one, and add the full NNLO cross section, as 


a 


NNLL 


{sH,Q^)=a^‘^^{sH,Q^) 


a 


’{sh, Q'^)\oiai) + 


a 


NNLO 


{sh, 


(17) 


and similarly for LL and NLL. 

With all the previous dehnitions and results, we are ready to perform the threshold resumma¬ 


tion up to NNLL. For more details of the resummation formalism, see for instance Ref. 26 


^'The lowest order cross section do starts at 0(a^ 
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Figure 1: The Higgs pair invariant mass distribution for Ecm = 14 TeV and the central scale /tq = Q, for 
the fixed order (left) and resummed (right) predictions. In the left (right) we show the LO (LL), NLO 
(NLL) and NNLO (NNLL) curves, with blue dotted, red dashed and black solid lines respectively. 




Figure 2: The Higgs pair invariant mass distribution for Ecm = 14 TeV and the central scale /tq = Q/2, 
for the hxed order (left) and resummed (right) predictions. The color coding is the same of Figure 

3 NNLL phenomenology 

We present in this section the phenomenological results. For the computation we take the Higgs 
mass to be Mh = 125 GeV. All the results are normalized by the exact LO top mass dependence, 
with Mt = 173.21 GeV. For the parton luminosities and strong coupling we use the MSTW2008 
sets, consistently at each perturbative order (i.e. LO PDFs and one-loop as evolution for LO 
and LL cross sections, etc.). The scale uncertainty was evaluated by varying independently the 
renormalization and factorization scales in the range /ro/2 < < 2/io with the constraint 

1/2 < flu/< 2, where po is the central scale. The analysis was performed for two choices of 
the central scale: = Q and po = Q/2, being Q the invariant mass of the Higgs pair system. 

The contributions from all the relevant partonic channels are always included in our numerical 
results. As described in the previous section, the threshold resummation only applies for the gg 
channel. With the corresponding matching we also account for the other partonic subprocesses at 
the corresponding hxed order accuracy. 

We start by showing the Higgs pair invariant mass distribution for a collider center of mass 
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Figure 3: The ii'-factors for the hxed order and resummed cross sections as a function of the Higgs pair 
invariant mass, for Ecm = 14 TeV. The left (right) panel shows the results for fio = Q = Ql‘^)- The 
color coding is the same of Figure 

energy Ecm = 14 TeV. In Figure we present the results corresponding to the central scale 
/To = Q, while in Figure [^the ones corresponding to /io = Q/2 are shown. For both figures, in the 
left plot we present the fixed order prediction (at LO, NLO and NNLO) while in the right one we 
show the resummed cross section (at LL, NLL and NNLL). 

In the hrst place we can observe that, with the exception of the /xq = Q/2 resummed distribu¬ 
tions, there is no overlap between the LO (LL) and NLO (NLL) bands, and it is only at second 
order that a sensible superposition of the bands occurs. We can also see from the plots that at 
every order the inclusion of the resummed contributions results in an increase of the cross section. 
Also, we can observe that the size of the uncertainty band at NNLL is always smaller than the 
corresponding NNLO one. This effect is more clear with the choice po = Q, for which also a better 
overlap between the NNLL and NLL bands is observed, with respect to the NNLO and NLO ones. 
The hxed order and resummed distributions have less differences for po = Q/2, as was already 
observed for single Higgs production, where the choice po = Mu/2 partially mimics some of the 
threshold resummation effects. Regarding the shape of the distributions, we observe very small 
differences after the resummation is performed. This is due to the fact that the relative size of 
the resummed contributions has a rather small dependence on the Higgs pair invariant mass. 

In Figure we present the iF-factors, dehned as the ratio between a given prediction and the 
LO one. For the denominator we hx /tr = fip = /iq. We observe, in more detail, the same features 
described above at the level of the cross section. In particular, it is visible that the resummed 
series has a better convergence than the fixed order one, exhibiting a larger overlap between the 

^For simplicity, we always label our resummed predictions as LL, NLL and NNLL. As explained before, these 
results include the matching to the fixed order cross section, so they should be interpreted as LL-l-LO, NLL-I-NLO 
and NNLL-I-NNLO respectively. 
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Figure 4: The ratio between the NNLL and the NNLO predictions as a function of the Higgs pair 
invariant mass, for the scales fi = Q (left) and /x = Q/2 (right). Results are shown for center of mass 
energies of 8 TeV (orange solid), 14 TeV (magenta dashed), 33 TeV (purple dot-dashed) and 100 TeV 
(black dotted). 
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Table 1: The total cross section and theoretical uncertainties for different center of mass energies, at 
NNLO and NNLL, for fj-o = Q and /xq = Q/2. PDF and PDF-|-as uncertainties correspond to the 
resummed predictions, and are estimated using the sets of MSTW2008 at 90% confidence level. 


first and second order bands. 

In Figure we show the ratio between the NNLL and the NNLO predictions, again as a 
function of the Higgs pair invariant mass, for different collider energies. The ratio shows an 
almost linear dependence on Q, increasing for higher invariant masses. Actually, this is expected 
because resummation contributions are enhanced when the process becomes closer to the partonic 
threshold. The same feature is reflected by the fact that the resummation contributions are 
relatively smaller for larger collider energies. We can also observe, as it was already clear from 
Figures and that the ratio between NNLL and NNLO is signihcantly smaller for the scale 
choice /xh = Rf = R = Q/2. At the total cross section level, for example, we hnd that the increase 
in the NNLL result with respect to the NNLO prediction is of 6.8% for Ecm = 14 TeV and fi = Q, 
while it drops down to 0.65% for fi = Q/2. 

We focus now on the theoretical uncertainty arising from the missing higher order contributions, 
which is estimated by the scale variation indicated above. In Table we present the total cross 
section predictions at NNLO and NNLL, together with the scale uncertainty. We can observe 
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Figure 5: The scale dependence of the total cross section at Ecm = 14 TeV, for the fixed order (upper) 
and resummed (lower) predictions. The color coding is the same of Figure 


that in all the cases the uncertainty of the resummed result is lower than the fixed order one. For 
instance, for Ecm = 14 TeV we hud that the total uncertainty at NNLO is 17% for /xq = Q, while 
it goes down to 11% at NNLL. This reduction is less important but still noticeable for /xq = Q/2, 
where it goes from 13% to 10% for the same center of mass energy. 


Another interesting feature to notice is the stability of the resummed cross sections, being 
almost independent of the choice of the central scale. The differences between the fio = Q and 
A^o = Q/2 results are below 2% in the central value, and present very similar uncertainty bands. 


In Table we also present the uncertainties coming from the strong coupling and parton flux 
determination for the resummed cross section. These were estimated using the MSTW2008 90% 
C.L. error PDF sets 31 , which are known to provide very close results to the PDF4LHC working 
group recommendation for the envelope prescription 32 . The results are very similar to the ones 
corresponding to the hxed order calculation (Tol. We can see that the PDF+as uncertainty is 


typically larger than the one arising from scale variation. 


In order to further illustrate the reduction of the uncertainty, we present in Figure the 
scale dependence of the total cross section, both for the hxed order and resummed predictions, 
at the different accuracy levels. The plots in the left correspond to varying simultaneously the 
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Figure 6: The total fixed order (FO) and resummed (RES) cross sections at leading (blue circle), next- 
to-leading (red square) and next-to-next-to-leading (black triangle) accuracy, for Ecm = 8 TeV (left) and 
14 TeV (right), for both central scales = Q and /tq = <5/2. The vertical solid lines indicate the scale 
uncertainty. The horizontal dotted lines indicate in each case the best prediction (NNLL). 


factorization and renormalization scales with fi^ = fip = the one in the center corresponds to 
the dependence on the renormalization scale for hxed /if = Q, while the plot in the right shows 
the factorization scale dependence, for hxed = Q. 


In all the plots we can observe that the inclusion of higher order corrections reduces the 
scale dependence. The contributions from resummation further reduce this dependence at NNLL, 
except for the fxp dependence at hxed fXR = Q. This last feature is also present for single Higgs 
production 26 , and suggests that the rather hat dependence on pi? at NNLO can be an accidental 
ehect. 


We summarize our results for the total cross section in Figure for Ecm = 8 TeV and 14 TeV. 
In this hgure we present the value of the total cross section at LO, NLO and NNLO along with 
the LL, NLL and NNLL predictions. The left panel of each plot shows the results corresponding 
to the central scale fxo = Q, while in the right one are those associated to po = Q/2. The vertical 
lines indicate the scale uncertainty of each result. 


The plots for both c.m. energies have similar features. As was pointed out before, for fxo = Q 
the corrections coming from threshold resummation are sizeable at every order, and the reduction 
of the total scale uncertainty at NNLL is notorious. For po = Q/2 the increase on the total cross 
section is much smaller, though one can still observe a reduction of the uncertainty. We can notice 
again the stability of the resummed prediction on the choice of the central scale. This is illustrated 
in Figure by the horizontal dotted lines, which indicates the NNLL result for each value of po- 
We can observe that the overlap between the two results is almost perfect, while in the case of 
the fixed order prediction there are much larger differences between them. 


Finally, we comment on the results obtained in Ref. 


33 


In that paper, a NLO-I-NNLL 
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prediction for Higgs pair production was presented, based on the soft-collinear effective theory. 
That calculation did not include the matching to the NNLO cross section nor the coefficient 
which were not available by the time of its publication. They found that the remaining total scale 
uncertainty was below 8%, a value that is slightly below our estimation. Also, their central value 
for the total cross section is below our result by about 2 — 3%. 


4 Conclusions 


We performe the soft-gluon threshold resummation up to NNLL accuracy, including consistently 
the matching to the NNLO cross section. We work in the large top-mass approximation, normal¬ 
izing our results by the exact LO dependence. 

We hnd that the resummation results in an increase of the total cross section of 6.8% for 
Ecm = 14 TeV and Hq = Q. The effect increases for lower center of mass energies and decreases 
for larger energies, as expected for the threshold contributions. The increase in the total cross 
section, for the same value of c.m. energy, goes down to 0.65% for /iq = Q/2. 


The scale uncertainty is also reduced with respect to the hxed order prediction, going from 
±8.5% to ±5.5%. The resummed prediction, including the corresponding uncertainty band, is 
found to be almost independent of the value chosen as the central scale, /iq = Q or /iq = Q/2. 
Given the similarity between the results, we can select in principle any of them as our hnal 
recommendation. For the sake of dehniteness, we recommend the usual setting = Q for the 
NNLL prediction, which for single Higgs production provides a result compatible with the recently 


computed N^LO cross section 34 


35 through the computation of 


36 via a reweighting technique that allows 


The hnite top-mass effects were analyzed at NLO in Ref. 
sub leading terms in the l/M^ expansion, and in Ref. 
to exactly include the one- and two-loop amplitudes. Based on those studies, we can estimate the 
uncertainties coming from the use of the effective theory to be ~ ±10% for the total cross section. 
The uncertainty coming from the missing higher orders of the QCD perturbative expansion is now 
dehnitely below that value, and also below the strong coupling and parton flux determination 
uncertainties, which is about ±7.5% for 14 TeV. 
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